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省 设 激 励 下 直升机 振动 咽 应 预计 及 实测 对 比 


王 形 ;于 文 湖 ;可 蝴 


(中 国 直 升 机 设计 人 研究 所 ,333001 景德镇 ) 


摘 要 : 旋 曙 激 励 是 直升机 振动 的 主要 源头 ,影响 直升机 驾驶 员 操 纵 性 、 机 上 人 员 和 舒适 性 及 结构 疲 
劳 寿命 。 振 动 问 题 一 般 在 试飞 过 程 中 暴露 出 来 ,问题 处 理 和 解决 往往 需要 大 的 周期 和 代价 ,影响 直 
升 机 研制 。 直 升 机 振动 响应 预计 是 解决 振动 问题 的 主要 手段 之 一 ,可 以 在 早期 发 现 振动 问题 ,方便 
在 型 号 研制 早期 开展 优化 设计 或 解决 预案 。 本 研究 基于 机 体 动 力学 仿真 模型 ,以 固定 坐标 系 下 的 
桨 虎 中 心 振动 载荷 作为 输入 ,得 到 机 体 关 键 位 置 的 振动 响应 ,形成 旋翼 激励 下 直升机 振动 响应 预测 
的 方法 。 通 过 菜 直 升 机 在 尾 交 激励 下 关键 位 置 响应 预计 与 实测 数据 对 比分 析 , 表 明 该 方法 得 到 的 
振动 响应 随 速 度 的 变化 规律 与 实测 数据 一 致 ,精度 满足 工程 需要 。 该 方法 在 新 机 振动 响应 分 析 、 振 
动 问题 排查 及 振动 控制 方案 制定 等 工作 中 具有 较 高 的 工程 价值 。 
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Prediction and comparing with measurement of helicopter 
vibration response excited by rotors 


WANG Ying, WANG Wentao, XI Juan 


(China Helicopter Research and Development Institute ,333001 Jingdezhen, China) 


Abstract :Rotor exitation is the main source of helicopter vibration ,which effects pilot manipulation ,com- 
fortability as well as structural fatigue. Vibration problems are generally exposed during test flights, and 
problem handling and solution often require large cycles and costs ,that affecting helicopter development. 
The prediction of helicopter vibration response is a major way to deal with vibration problems ,which can 
be found in the early stage ,convenient to carry out optimization design or solution plan in the early devel- 
opment stage of the helicopter. In this paper,a method to predicting vibratory response excited by rotors 
was presented , which takes vibratory loads of rotor center in nonrotating coordinate system as input , and 
based on the structural dynamic model. The method was proved to be effective ,whose results had coordina- 
ted variation patterns as measured data, as well as satisfying accuracy of the engineering needs , through 
comparing analysis of predicted vibration response and measured ones excited by tail-rotor of the important 
positions of a certain helicopter. The method has high engineering value in vibration response analysis ,Vi- 


bration problem investigation and vibration control scheme formulation. 
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直升机 振动 问题 是 直升机 研制 过 程 中 一 直 存在 
的 问题 ,影响 机 上 人 员 舒 适 性 .工效 性 ,机 上 设备 的 
功能 与 性 能 ,以 及 结构 的 疲劳 强度 。 如 何 解决 与 处 
理 好 振动 是 直升机 研制 成 败 的 关键 因素 之 一 。 

旋翼 系统 旋转 产生 的 周期 激 振 力 是 直升机 振动 
的 最 主要 源头 ,也 是 出 现 振动 问题 的 最 主要 原因 。 
常规 的 单 旋翼 带 尾 桨 构 型 直升机 , 尾 桨 作为 关键 部 
件 之 一 ,为 直升机 提供 反 扭矩 ,实现 航向 操纵 及 航向 
稳定 性 "““” 。 尾 桨 作为 直升机 高 速 旋转 部 件 , 是 直 升 
机 振动 的 主要 振 源 之 一 。 尾 桨 产生 的 振动 对 尾部 结 
构 及 部 件 的 工作 环境 产生 较 大 的 影响 ,国内 近年 来 
发 生 了 多 起 由 于 尾 桨 振动 过 大 导致 尾部 结构 裂纹 、 
摇 避 断裂. 脚 路 打脚 . 尾 桨 操纵 失效 等 现象 ,甚至 出 
现 尾 桨 断裂 导致 直升机 除 毁 的 事故 。 实 现 直升机 尾 
部 低 振 动 水 平 的 设计 目标 显得 尤为 重要 。 
加 以往 型 号 一 般 只 在 直升机 试飞 后 对 机 上 振动 水 
平 开展 实测 ,并 逐步 暴露 振动 问题 。 而 振动 问题 的 
代理 一 般 需要 花费 较 大 的 周期 和 代价 ,影响 型 号 的 
研 恒 进度 。 另 一 方面 ,解决 振动 问题 或 实现 低 的 振 
惑 目标 ,往往 需要 开展 旋翼 系统 的 优化 设计 ,或 针对 
振 沽 响应 位 置 等 采取 振动 控制 措施 ,在 试飞 阶段 很 
难 更 改 设计 或 采取 振动 控制 措施 。 相 比 而 言 ,能 
快速 准确 预测 尾 桨 振动 响应 , 在 型 号 研制 初期 掌握 
尾部 的 振动 环境 ,制定 振动 控制 预案 对 保障 直升机 
的 Ck 行 安 全 意义 重大 "7]。 
志 本 研究 基于 机 体 动力 学 模型 ,利用 CAMRAD I 
软 循 ,选取 典型 飞行 状态 进行 计算 ,得 到 桨 载 中 心 
NQ(N 为 桨 叶 数 量 ,0 为 旋 避 转速 ) 的 振动 载荷 ,以 
桨 坑 中 心 振动 载荷 作为 输入 ,得 到 机 体 关键 部 位 的 
振动 响应 ,形成 旋翼 激励 下 直升机 振动 响应 预测 的 
方法 ,实现 在 型 号 设计 初期 阶段 预计 振动 响应 的 
目的 。 通 过 某 直 升 机 在 尾 桨 激励 下 的 关键 位 置 响应 
预计 与 实测 数据 对 比分 析 , 表明 该 方法 得 到 的 振动 
响应 精度 可 接受 , 随 速度 的 变化 规律 与 实测 数据 一 
致 ,在 新 机 振动 响应 分 析 ,振动 问题 排查 及 振动 控制 
预案 等 工作 中 具有 较 高 的 工程 价值 。 


1 动力 学 建 模 


直升机 机 体 处 的 振动 决定 于 旋 踊 激励 载 答 和 旋 
以 至 机 体位 置 的 传递 函数 ,而 准确 的 建立 机 体 结构 


动力 学 模型 ,是 得 到 传递 函数 的 关键 。 

采用 有 限 元 方法 建立 机 体 动 力学 分 析 模 型 ,该 
模型 包括 除 旋 恤 以 外 的 结构 ,按照 机 体 的 主 传 力 路 
线 及 重量 重心 进行 建立 ,真实 地 模拟 全 机 的 刚度 和 
质量 分 布 。 根 据 研 究 , 采 取 以 下 原则 能 够 得 到 更 准 
确 的 机 体 动力 学 模型 。 

1) 结构 模型 中 ,对 机 体 的 各 个 框 纵 梁 、 地 板 、 蒙 
皮 、 短 轰 平 尾 、 上 、 下 垂 尾 等 主要 厌 力 结构 以 梁 单 
元 \ 完 单元 来 模拟 ,进行 精细 化 建 模 。 

2) 对 于 发 动机 、\ 主 减速 右 两 大 动 部 件 , 以 及 相应 
的 连接 支撑 组 件 ,参照 实际 结构 建立 实体 模型 。 

3) 对 于 前 、 后 罗 驶 员 座 椅 、 尾 传动 轴 、 旋 辟 、 尾 
浆 自动 倾斜 锅 、 桨 发 油箱、 外 挂 武器 .前 后 驾驶 员 、 
起 落架 以 及 机 载 系统 设备、 装置 等 ,采用 集中 质量 
单元 进行 模拟 ,将 这 些 集中 质量 根据 实际 的 安装 位 
置 情况 ,通过 刚性 连接 单元 或 加 权 平 均 约束 单元 连 
接 到 各 自 的 实际 安装 节点 上 。 

4) 机 号 结 构 中 的 连接 (螺栓 连接 、 锦 钉 连 接 等 ) 
大 部 分 采用 了 常用 的 节点 缝合 ,而 部 分 结构 之 间 的 连 
接 ,如 短 豆 与 机 且 . 平 尾 与 尾 梁 , 挂 架 与 短 细 武吉 与 挂 
架 等 , 则 根据 相应 的 传 力 形式 与 连接 性 质 ,分 别 选 择 采 
用 刚性 连接 单元 和 加 权 平 均 约束 单元 进行 模拟 。 


2 旋翼 动 载 计算 


2.1 载荷 的 来 源 


直升机 各 片 桨 叶 的 振动 载荷 通过 浆 载 传 向 机 
体 , 引 起 机 体 振动 。 根 据 旋 翼 动 力学 相关 理论 , 桨 载 
振动 载荷 的 来 源 主 要 有 浆 叶 根部 挥舞 剪 力 挥舞 弯 
和 抢 、 摆 振 剪 力 、 摆 振 弯 矩 及 扭矩 ,各 片 桨 叶 的 载荷 在 
桨 载 中 心 合成 ,转换 到 固定 坐标 系 下 对 应 浆 载 六 力 
素 , 即 拉力 水平 力 \ 侧 向 力 、 俯 仰 力矩 、 滚 转 力矩 . 捏 
和 矩 , 桨 根 载荷 (旋转 坐标 系 ) 与 桨 载 载荷 ( 固定 坐标 
系 ) 的 定义 分 别 如 图 1、 图 2 所 示 。 由 于 浆 载 的 “过 
滤 ” 效 应 , 仅 保留 了 kNQ(N 为 桨 叶 数 量 ,上 为 正 整 
数 ) 的 激 振 力 , 桨 叶 根 部 载荷 (旋转 坐标 系 ) 与 桨 载 
中 心 振动 载荷 (固定 坐标 系 ) 之 间 转换 见 式 (1) ~ 
(6) ,直升机 浆 叶 根部 高 阶 ((NW+1)O、NO) 载 荷 ， 
经 桨 坑 合 成 传递 到 机 体 结构 上 ,形成 机 体 NO 振动 
的 主要 激励 源 , 桨 根 载荷 与 桨 载 中 心 振 动 载荷 对 应 
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关系 总 结 如 图 2 所 示 。 其 中 : 垂 向 力 ( 拉 力 ) 了 来 源 


于 桨 叶 根 部 挥舞 方向 ( 垂 向 ) NO 的 剪 力 ; 水 平 力 H、 
侧 向 力 了 来 源 于 浆 叶 根部 旋转 面 内 ( 摆 振 方向 和 展 
向 )(N+t1)Q 的 前 力 ;俯仰 力矩 W 滚 转 力矩 ML 来 
源 于 奖 叶 根 部 (N +1)0Q 挥舞 弯 算 ;扭矩 0 来 源 于 桨 
叶 根 部 NQ 氛 振 弯 和 矩 。 


T= > NS ve (1) 
p=-% 
H= A 1S +iS, Jew™ 
SS 2 TpN-1 TpN+1 xpN-1 XpN+1 
(2) 
三 SS VN ipNy 
= 2 2 (a iS, N-1 二 LS XpN-1 XpN+1 
(3) 
有 oo 
0= 》 NN,, em (4) 
P=-% 


= N ip! 
M, > 7\ | Ne Nps ) e" (5) 
P=-% 
N 和 IPN 
二 二 之 ， 本 -TV +i Nps)e ” (6) 


: 了 为 桨 开 垂 向 力 (拉力 ) ; 为 桨 红 水 平 力 ; 了 
J 尝 载 侧 向 力 ; @ 为 桨 坑 扭 矩 ; M, 为 桨 载 深 转 力矩 ; 
1 为 桨 坑 俯 仰 力矩 ; N 为 奖 叶 数量 ; 少 为 第 m 片 桨 
只 的 方位 角 ; 5, 为 奖 叶 根部 挥舞 方向 剪 力 ; 5, 为 浆 
叶 浪 部 径 向 剪 力 ; 5, 为 桨 叶 根部 摆 振 方向 前 力 ; N， 
为 奖 叶 根部 摆 振 力矩 ; N; 为 桨 叶 根 部 挥舞 力矩 。 


AN 
CY > 
= 
J 
S. 


图 1 桨 根 载荷 定义 (旋转 坐标 系 ) 


Fig.1 Rotor forces and moments acting 


on the hub (rotating frame) 


图 2 桨 载 振 动 载荷 定义 (固定 坐标 系 ) 


Fig.2 Rotor forces and moments acting 


on the hub (nonrotating frame) 


(hinaX i /全 人 FE 甘 日 干 
工 影 ， 等 : 旋 辟 激 励 下 直升机 振动 响应 预计 及 实测 对 此 ninaXiv 人 富 1 人 ET 255 


表 1 桨 根 载荷 与 桨 坑 振 动 载荷 对 应 关系 


Relationships between root loads and hub vibration loads 


Tab. 1 
桨 根 载荷 (旋转 坐标 系 ) 


桨 载 振动 载荷 (固定 坐标 系 )pNQ 


7 5S.; pNQ 
0 Ni; PNQ 
HY SS; (pN+1)0 
M, M. Ne;(pN+1)0 


2.2 ” 奖 载 振动 载荷 计算 


CAMRAD 本 是 一 款 旋 浑 飞行 器 综合 分 析 软 件 ， 
包含 多 体 动力 学 , 非 线 性 有 限 元 结构 动力 学 和 旋翼 
空气 动力 学 等 分 析 模 型 。CAMRAD 工 可 用 于 多 个 
旋翼 飞行 器 构 型 (如 孤立 旋 避 . 单 旋 辟 带 尾 桨 直 升 
机 、 倾 转 旋 避 机 、 共 轴 直 升 机 、 多 旋翼 直升机 ) ,可 进 
行 性 能 、 载 荷 、 振 动 、 动 响应 、 稳 定性 多 种 计算 分 
析 " ,一 般 的 计算 框架 如 图 3 所 示 。 


全 总 和 1 r pt 区 1 

! 气动 力 模型 ! ! ”部 件 模 型 ! 

a 
入 流 模型 旋翼 


1 

' 

l 天 
站 机 身 Ts ; 

1 

1 


图 3 CAMRAD 开 综合 分 析 框 架 
Fig.3 Framework of CAMRAD JI comprehensive analysis 

CAMRAD 本 气动 模型 主要 基于 升力 线 理 论 和 
动量 理论 ,主要 包括 定常 人 流 模型 .动力 人 流 模 型 和 
自由 尾 迹 模型 等 ,用 于 不 同 用 途 下 的 诱导 速度 .动态 
失速 及 尾 迹 描述 。 结 构 模 型 主要 包括 刚体 梁 模 
型 和 弹性 梁 模 型 ,用 于 不 同 用 途 下 的 结构 特性 描 
述 '"”。 定 常 入 流 模 型 及 刚体 梁 模 型 可 用 于 计算 性 
能 、 功 率 需 求 等 ,模型 简单 ,收敛 速度 快 ;复杂 的 气动 
模型 及 弹性 深 模 型 可 用 于 计算 旋翼 载荷 .振动 及 噪 
声 等 ,模型 复杂 ,不 易 收 敛 。 

由 于 浆 载 振动 载荷 (固定 坐标 系 下 ) 是 由 桨 叶 
高 阶 ((CY+1)C2、NO) 谐 波 载 荷 琶 加 而 成 (转换 关系 
见 式 (1) ~ (6))。 因 此 借助 CAMRAD 开 计 算 浆 载 
振动 载荷 时 , 需 采 用 自由 尾 迹 模型 以 及 弹性 梁 模 型 ， 
其 中 自由 尾 迹 模型 考虑 了 浆 叶 尾 迹 引 起 的 非 均 匀 诱 
导 速 度 的 影响 ,弹性 梁 模 型 考虑 了 桨 叶 高 阶 模 态 在 
气动 激励 下 的 响应 。 例 如 ,不 考虑 桨 叶 尾 迹 影 响 时 ， 
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某 型 机 MO 垂 向 力 小 于 50 N , 而 采用 自由 尾 迹 计算 
模型 时 ,NO 垂 向 力 约 为 3 000 N ,更 接近 工程 实际 。 


3 振动 响应 预计 


旋 喜 激励 下 的 振动 响应 是 机 体 结 构 振 动 的 主要 
部 分 ,其 一 般 的 传递 路 径 为 : 桨 叶 一 桨 载 一 减速 絮 一 
机 体 结构 一 机 体 各 部 件 ,其 量 值 决 定 于 桨 载 振动 载 
荷 及 桨 载 中 心 到 各 位 置 处 的 传递 函数 。 浆 载 中 心 振 
动 载 位 激励 下 关键 位 置 的 振动 响应 计算 关系 式 见 式 
(7)。 借 助 全 机 动力 学 模型 ,在 桨 融 中 心 处 加 载 桨 


坊 振 动 载荷 ,获取 关注 位 置 的 振动 响应 。 振 劲 响应 
计算 流程 图 见 图 4。 


H,. 
xX 
H Qy ( 7 ) 
万 和 
Oz 0 


K 轩 ; u,、u，、u, 为 关注 位 置 3 个 方向 的 振动 响应 ; 
且 夺 对 应 桨 载 中 心 某 个 方向 的 载荷 对 关注 位 置 菜 个 
] 的 传递 阴 数 ;TX 等 对 应 浆 融 中 心 6 个 方向 的 


y 


J 载 和 从。 


弹性 梁 
模型 


桨 壳 振 动 载荷 


动力 学 模型 
振动 响应 


图 4 振动 响应 计算 流程 图 


Fig.4 Flowchart of vibration response calculation 


4 实例 验证 


以 某 直升机 尾 浆 安装 处 振动 响应 计算 为 例 , 开 
展 了 在 尾 桨 激励 下 的 关键 位 置 啊 应 预计 与 实测 数据 
对 比分 析 。 

依据 直升机 三 维 结构 数 模 ,按照 第 1 章 所 述 原 
则 ,建立 某 直 升 机 全 机 动力 学 模型 ,模型 中 的 单元 统 
计 见 表 2 , 主要 模 态 计算 结果 与 试验 结果 对 比 见 
表 3。 

利用 CAMRAD 本 软件 ,选取 典型 飞行 状态 进行 
计算 , 桨 载 中 心 NO 的 振动 载荷 如 表 4 所 示 。 将 表 4 
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1 | 
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的 尾 桨 浆 纹 中 心动 载荷 变化 趋势 绘制 如 图 5 所 示 。 

得 到 尾 桨 安装 处 航线 (X)、 侧 向 (2 切 及 垂 向 (Z) 的 振 

动 响应 仿真 数据 如 表 5 所 示 。 将 表 5 的 尾 桨 安装 处 

的 振动 响应 仿真 数据 变化 趋势 绘制 如 图 6 所 示 。 
表 2 动力 学 模型 单元 统计 


Tab.2 Elements statistic of dynamic model 


类 别 ”节点 杆 单 元 ” 梁 单 元 ”这 单元 ”刚体 元 质量 单元 


数量 ”46 296 20 10986 44942 37 210 


表 3 模 态 计算 结果 对 比 


Tab.3 Comparison of calculation and measurement 


计算 结果 / ”试验 结果 / 相对 误差 / 下 

下 和 2 模 态 名 称 

1 7.8 7.9 -1.3 向 一 阶 
2 8.7 8.8 -1.1 垂 向 一 阶 
3 14.9 15.0 -0.7 出 向 二 阶 
4 14.6 14.4 1.4 垂 向 二 阶 
5 25.3 24.3 4.1 岗 向 三 阶 
6 26.0 26.1 -0.4 垂 向 三 阶 


表 4 尾 桨 席 振动 载荷 
Tab.4 Vibration loads of tail rotor hub 
速度 / F/ FF/ FF/ M.,/ M,/ M./ 


(km*:h-!) N N N (N-:m) (N:m) (N:m) 


80 325.0 222.4 354.1 326.6 214.2 350.6 


110 281.5 235.7 302.4 351.2 256.4 381.6 


140 215.2 249.2 243.9 340.2 262.2 376.4 


180 148.0 293.8 160.6 315.0 219.2 349.5 
200 114.2 258.0 131.7 294.8 223.3 321.6 
230 188.3 268.2 210.6 386.6 431.6 405.7 


250 261.0 313.0 274.5 460.5 590.4 481.3 


扭矩/(N ，m) 


50 100 150 200 250 300 
速度 /dm * hh ") 


Fig.$5 Vibration loads of tail rotor hub 
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表 5 尾 桨 安装 处 振动 响应 仿真 数据 规律 一 致 , 幅 值 平均 误差 约 为 20% ,结果 相对 可 信 ， 
Tab.5 Vibration response of tail rotor 可 满足 工程 需要 。 
速度 /lm'hiD ee au 3) 计 算 结果 与 实测 振动 水 平 之 间 的 误差 ,主要 

g0 5 下 Tt 来 源 于 桨 载 振动 载荷 计算 的 复杂 性 ,动力 学 模型 的 

二 | a 简化 以 及 响应 的 非 线性 。 

140 1.42 0.75 0.93 

180 1.34 0.96 0.76 

200 1.25 1.26 0.72 

230 3;15 1.93 1.95 

250 3.70 2:21 2.44 


50 100 150 200 250 300 
速度 /mh ) 


5 多 8 尾 桨 安装 处 振动 响应 仿真 与 实测 数据 对 比 -X 
玉 Fig.8 Comparison of measurement vibration response 
< 二 号 with simulation of tail rotor-X 
©O 25 
©O | 
©O 20 

50 100 150 200 250 300 
es 速度 /(km i") Lo 15 人 

疾 

3 图 6 尾 桨 安装 处 振动 响应 和 加 
© Fig.6 Vibration response of tail rotor 二 
EN 某 型 机 在 科研 试飞 阶段 ,开展 了 不 同 飞行 状态 2 


下 葛 振 动 水 平 测试 ,采用 快速 依 里 叶 变 换 处 理 得 到 ‘50 10 150 20 250 300 
ys 这 是 a i 速度 /(km。h 
尾 涵 安 装 处 NO 振动 水 平 ,实测 振动 水 平 随 飞 行 束 CO 
度 恶化 如 图 7 所 示 , 各 个 方向 计算 结果 与 实测 数据 加 9 尾 桨 安装 处 振动 响应 仿真 与 实测 数据 对 比 - 
对 入 如 图 8 ~10 所 示 Fig.9 Comparison of measurement vibration response 
a with simulation of tail rotor-Y 
3 
© 
3.0 
2.5 这 
-0 2.0 
四 2 Lo 
二 1.5 加 
1.0 
下 
65 
0.0 
50 100 150 200 250 300 
速度 /(km .1 50 100 150 200 250 300 
速度 /dm .h ) 
将 二 着 均 等 溃 涩 i > 局 > 四 | 滥 | 
了 图 10 尾 桨 安装 处 振动 响应 仿真 与 实测 数据 对 比 -Z 
Fig.7 Measured vibration response of tail rotor Fig. 10 ”Comparison of measurement vibration response 
通过 响应 的 计算 结果 及 与 实测 结果 对 比分 析 ， with simulation of tail rotor-Z 
主要 规律 如 下 。 


1) 尾 桨 安装 处 振动 响应 随 速 度 增 加 先 减 小 后 增 SS 结 论 
大 ,速度 大 于 200 km/h 后 各 方向 幅 值 显著 增加 。 
2) 振 动 响应 计算 结果 与 实测 振动 数 随 速 度 变化 本 研究 基于 机 体 动力 学 模型 ,以 固定 坐标 系 下 
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的 桨 载 振动 载 谷 作为 输入 ,进行 啊 应 预测 分 析 ,得 到 
机 体 关 键 位 置 的 振动 响应 ,形成 旋 泪 激励 下 直升机 
振动 响应 预测 的 方法 ,实现 了 在 型 号 研制 初期 阶段 
预测 关键 位 置 振动 响应 的 目的 ,形成 主要 结论 如 下 。 

1) 通 过 某 直升机 在 尾 桨 激励 下 关键 位 置 响应 预 
计 与 实测 数据 对 比分 析 , 表 明 该 方法 得 到 的 振动 响 
应 随 速度 的 变化 规律 与 实测 数据 一 致 , 幅 值 精度 满 
足 工 程 需要 。 

2) 该 方法 在 新 机 振动 响应 分 析 、 振 动 问题 排查 
及 振动 控制 方案 制定 等 工作 中 具有 很 高 的 工程 
价值 。 
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